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( The fascinating electronic properties of the family of layered organic molecular crystals k-(BEDT- 

TTF) 2 X where X is an anion (e.g., X=I 3 , Cu[N(CN) 2 ]Br, Cu(SCN) 2 ) are reviewed. These materials 
are particularly interesting because of similarities to the high-T c cuprate superconductors including 
unconventional metallic properties and competition between antiferromagnetism and superconduc- 
tivity. The temperature dependence of electrical transport, optical, and nuclear magnetic resonance 
properties deviate significantly from those of a conventional metal. In particular, there appears to 
be an effective Fermi energy of the order of 100 K which is an order of magnitude smaller than pre- 
dicted by band structure calculations. The results of quantum chemistry calculations suggest that a 
minimal theoretical model that can describe these materials is a Hubbard model on an anisotropic 
triangular lattice with one electron per site. Appropriate parameter values for the model imply that 
the electronic correlations are strong, significant magnetic frustration is present, and the system is 
close to a metal-insulator transition. Insight into the physics of this model can be obtained from 
recent studies of the Hubbard model using a dynamical mean-field approximation. They are con- 
sistent with a low effective Fermi energy and the unconventional temperature dependence of many 
of the properties of the metallic phase. Future directions are suggested for both theoretical and 
experimental studies. 

a 

The discovery of heavy fermion metals, high-temperature superconductivity in copper oxides, and colossal magne- 
i O | , toresistance in manganates has stimulated extensive experimental studies of these materials and extensive theoretical 
studies of strongly correlated electron models £1 A key feature of the copper oxides (cuprates) is their layered or 
| quasi-two-dimensional structure which enhances the effects of the -interactions between the electrons and leads to 
■ metallic properties that are quite distinct from conventional metals.Hu Organic conductors based on the BEDT-TTF 
[bis-(ethylenedithia-tetrathiafulvalene)] molecule have been widely studied because prior to the discovery of fullerene 
superconductors they had the highest transition temperature (T c ~ 10 K) of an organic superconductor. The fam- 
ily k— (BEDT-TTF)2X whes&X is an anion (X=Cu[N(CN)2]Br, for example) are particularly interesting because of 
similarities to the cuprates an This article gives a brief review of these materials focussing on how the properties of 
qq ■ the metallic phase are quite different from those of a conventional metal. The, main purpose of this article is to argue 
that this is because these materials are strongly correlated electron systemsQ and that a model introduced by Kino 
"j—t ■ and FukuyamaQ may be the simplest possible theoretical model to describe them. Hopefully the article will stimulate 
more experimental and theoretical work on these exciting materials- Readers interested in a broader review of organic 
superconductors should consult recent monographs on the subject.ETtlj 
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I. EXPERIMENTAL PROPERTIES 

Organic molecular crystals with the formula (BEDT-TTF^X consist of conducting layers of BEDT-TTF molecules 
sandwiched between insulating layers of anions X. The layered structure leads to highly anisotropic electronic proper- 
ties. BEjDT-TTF is a large planar molecule and the different possible packing patterns are denoted by different Greek 
lettersO. The basic unit of the packing pattern in the k phase is a "dimer" consisting of two BEDT-TTF molecules 
stacked on top of one another. Each dimer has one electron less than a full electronic shell because of charge transfer 
to the anions. 

The family k— (BEDT-TTF^X has a particularly rich phase diagram as a function of pressure, temperature, and 
anion X=Cu[N(CN) 2 ]Cl, Cu[N(CN) 2 ]Br, and Cu(SCN) 2 , as shown in Fig. 1. A number of features of this diagram 
should be noted, (i) Antiferromagnetic and superconducting phases occur next to one another, (ii) At low pressures 
the metallic phase has properties that are quite distinct from conventional metals, (iii) The phase diagram is quite 
similar to that of the cuprates if pressure is replaced with doping. The properties of the different phases are now 
reviewed. 
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A. Metal-insulator transition 



In X=Cu[N(CN)2]Cl under moderate pressureS|420-30 Mpa) as the temperature is lowered there is a first-order 
transition from an insulating into a metallic phaselij This transition also occurs at constant tempeijatures between 10 
and 30 K as the pressure is increased. At ambient pressure X= Cu(CN) 3 is a paramagnetic insulator £3 For pressures of 
0.35 to 0.4 GPa there is a first order transition as the temperature is lowered. (For 0.35 GPa the resistivity decreases by 
five orders of magnitude at the transition) . The temperature at which the metal insulator transition occurs increases 
with pressure. It will later be argued that this metal-insulator transition is a Mott-Hubbard transition driven by 
electron-electron interactions. However, it needs to be checked whether there is any structural change associated with 
this transition. V2O3 undergoes a metal-insulator transition. Since there is a jump in one of the lattice constants at 
the transition its origin could be structural rather than electronic. 



B. Insulating antiferromagnetic phase 

Evidence of a&tiferromagnetic ordering below 26 K is provided in X=Cu[N(CN) 2 ]Clll and deuterated 
X=Cu[N(CN) 2 ]BiO by the observed splitting of proto n nmr lines. The magnetic moment estimated from the magni- 
tude of the splitting is (0.4 — 1.0) (13 per dimer. This large moment suggests that there are strong electron-electron 
interactions in these materials. Significant anisotropy is seen in the magnetic susceptibility below.-26 K. Below 22 K 
there is weak ferromagnetic hysteresis with a saturation moment of about 10~~ 3 hb pepjjormula unit.E3 This small ferrp-. 
magnetic moment could be due to a slight canting of the antiferromagnetic moments.E-3 Antiferromagnetic resonanceta 
results are consistent with an easy plane antiferromagnet with Dzyaloshinskii-Moriya interaction. The antiferromag- 
netic phase becomes unstable under moderate pressures of about 300. bar. It can be stabilised and enhanced by a 
magnetic field perpendicular to the layers of the order of a few tesla.Ea 

The observed phase diagram at the boundary of the superconducting and insulating phases is actually more com- 
plex than shown in Fig. 1. For example, coexistence of superconductivity and antiferromagnetic phases has been 
observed. E3E3Ej As the temperature is lowered a superconductor to insulator transition is observed. 



C. Metallic phase 



Many of the properties of the metallic phase have a temperature dependence that is quite distinct from that of 
conventional metals. Yet at low temperatures (less than about 30 K) some properties are similar to those of a 
conventional metal, but with a Fermi energy of the order of 100 K. This is almost an order of magnitude smaller 
than predicted by band structure calculations. Application of pressures of the order of 10 kbar restores conventional 
metallic properties over the full temperature range. The details of different properties are now reviewed. 

Optical conductivity. InfraredElrEj and microwaveEj measurements of the frequency dependent conductivity <j(uj) 
of X=Cu(SCN)2, X=Cu[N(CN)2]Br, and X=i3L3 deviate from the Drude behavior found in conventional metals. At 
room temperature a(u) is dominated by a broad peak around 300 or 400 meV (depending on the polarization and 
anion X) with a width of about 150 meV. Even down to 50 K no "Drude-like" peak at zero frequency is present (see 
Fig. [2|). At 25 K the high energy peak decreases slightly in temperature and a "Dcude-like" peak but can only be fit 
to a Drude form if the scattering rate and effective mass are frequency dependent .Ell 

Nuclear magnetic resonance. In a conventional metaL the Knight shift K is proportional to the density of states at 
the Fermi energy and is independent of temperature. Ej In contrast, in X=Cu[N(CN)2]Br the Knight shift decreases 
significantly below about 50 K, suggesting a suppression of the density of states or Jiipseudogap" near the Fermi 
energy.ESEH In conventional metals the NMR relaxation rate 1/Ti obeys Korringa's law:E£l l/(TiTK 2 ) is independent 
of temperature and a universal number. The relaxation rate is five to ten times larger than expected from Korringa's 
law and is strongly temperature dependent: L/.C24X) has a peak in the range 10 to 50 K, depending on how close 



the system is to the metal-insulatc 



MsitionEa€3Ej Above 50 K the relaxation rate is similar for X=Cu[N(CN) 2 ]Cl, 
Cu[N(CN)2]Br and X=Cu(NCS)2E3't£l. The abepee behavior has some similarities to the temperature dependence of 
the nmr properties of the underdoped cuprates.Ej As the pressure is increased to four kilobars the NMR properties 
become more like those of a conventional metal: the Knight shift and \/{TiT) become independent of temperature.EZI 
Temperature dependence of the resistance. Unlike a conventional metal the resistance does not monotonically 
increase with temperature. It is a maximum at T max where T™- ~ 100 K for X=Cu(NCS)2- With increasing 
pressure T max increases and the peak disappears at high pressures 113 From abottt 30 K down to the superconducting 
transition temperature the resistance decreases quadratically with temperatureEll 
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p(T) = po + AT 2 (1) 

Such a temperature dependence is characteristic of metals in which the dominant scattering mechanism is the in- 
teractions of the electrons with qae-another and is observed in transition metals and heavy fermions. In those 
systems the Kadowaki- Woods rulet3H relates the coefficient A in ([[]) to the linear coefficient for the specific heat, 7: 
A/7 2 =constant. The constant is 4.0 x I0~ 13 51cm (mol/mJ) 2 for transition metals, and-,1.0 x 10 -11 fkm (mol/mJ) 2 
for heavy fermions and for transition metal oxidesHiear the Mott-Hubbard transition! 3 -! In the organics the ratio 
is 5 to 200 times larger than predicted by this IotHeJ. It is has been suggested that this discrepancy implies that 
the daminant scattering mechanism is phonons.BEj In some strong coupling superconductorstiH and dirty metallic 
filmsE3 a T 2 resistivity is observed. If interference between electron-phonon scattering and disorder is taken into 
account theory predicts that A ~ ap(Q)/6 2 D where a ~ 0.01 — 0.1, and p(0) is the residual resistivity due to disorder. 
Comparing to the data on X=Cu(SCN)2 in Ref. 5lJ the value of A predicted by this expression is several orders of 
magnitude smaller than observed. I— ■ 

Mean-free path and Mott's minimum conductivity^ If the electronic mean- free path is less than a lattice constant 
than the idea of electronic transport by electrons with well-defined wavevector is no longer meaning ful.El Then the 
material should act as an insulator with a resistance that drops with increasing temperature; to the contrary it remains 
a metal with an increasing resistance. If the mean-free path is of the order of a lattice constant then for a quasi-two 
dimensional Fermi surface the intralayer conductivity will be approximately, 

0~min ^ , \") 
h C 

where c is the interlayer spacing and B is a constant of order one. For c ~ IOA this gives o~ m i n ~ 10 3 (f2cm) — 1 . The 
intxalayer conductivity of X=Cu(NCS)2 is comparable to this at 30 K and is three orders of magnitude smaller at 100 
K.El However, from 30 to 100 K, the resistance increases with temperature, characteristic of a metal. Hence, these 
materials can be classified as "bad metals. "Hi 

Hall resistance. In a conventional metal the Hall resistance is a measure of the pumber of charge carriers and is. 
independent of temperature above about 0.3#d, where 9d is the Debye temperature. c3 (Specific heat measurementsEl 
suggest Op ~ 200 K in the K salts). The Hall resistance of X=Cu(SCN) 2 is weakly temperature dependent above 80 K 
but increases by a factor of three as the temperature is decreased from 80 K down to the superconducting transition 
temperature (10 K)H Hall measurements on X=Cu[N(CN) 2 ]Br and X=Cu[N(CN) 2 ]Cl at ambient pressures also 
found a strong temperature dependence but were complicated by long-term sample resistance relaxation. Recent 
measurements on the X=Cu[N(CN)2]Cl salt under pressure also found a strong temperature dependences. However, 
the ratio of the longitudinal resistance R xx to the Hall resistance Rh has a simple quadratic temperature dependence 
up to 100 K: 

cot6 H = ^ = aT 2 + C. (3) 

9h is known as the Hall angle. Similar behavior is found in the cuprate superconductorsdi. It has been suggested .that 
in the cuprates that there are two distinct scattering times and that these have a different temperature dependence e3'c3 
The Hall resistance involves the product of the two scattering rates. This hypothesis needs to be tested in the organics 
by checking if R 2 x /Rh is independent of temperature. 

Thermoelectric power. In a conventional metal the thermoelectric power is given by 

^2 J.2 rp 

(4) 
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where fc^ is Boltzmann's constant, e is the electronic charge, and Ep is the Fermi energy. The Fermi energy estimated 
from comparing this expression to measurementsEJ on X=Cu[N(CN)2]Br is several times smaller than predicted by 
Hiickel band structure calculations. Furthermore, the thermopower is not linear in temperature above 30 K and 
has a peak around 100 K (as does the resistivity) and is similar in magnitude and temperature dependence to the 
cupratesE2l. 

Specific heat. The low-temperature specific heat of X=Cu[N(CN)2]Br and X=Cu(SCN)2 have been measureoS 
at ambient pressure and under a magnetic field of 10 tesla (to destroy the superconductivity). The specific heat 
coefficients 7 are 22 ± 3 mJ K -2 mol _1 |-and 25 ± 3 mJ K~ 2 mol -1 , respectively, about 2.5 times larger than predicted 
by Hiickel band structure calculations.E 2 ! 

Magnetic susceptibility. In a conventional metal this is independent of energy for temperatures much less than than 
the Fermi energy. Fpf, the k salts the susceptibility \ is almost independent of temperature above 50 K but decreases 



slightly below 50 K 
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Wilson ratio. The dimensionless quantity R = 47r 2 fc^x(0)/(3(ffAis) 2 7) where g ~ 2 is the gyro-magnetic ratio and 
fiB is the Bohr magneton. For a non-mteracting Fermi gas R — 1 and for the Kondo model R = 2. For heavy fermion 
metals it is in the range one to three.^ For both X=Cu[N(CN) 2 ]Br and X=Cu(SCN) 2 it is 1.5 ± 0.2. 

Magneto-oscillations. The temperature dependence of the amplitude of Shubnikov - de Haas and de Haas - van 
Alphcn oscillations can be used to deduce the effective mass, m*, of the charge carriers associated with various 
orbits on the Fermi, surface. □ Values, obtained £p£ the fj orbit include m*/m e = 5.4 ± 0.1,7.1 ± 0.5,3.9 ± 0.1 for 
X=Cu[N(CN) 2 ]Br,pLX|=Cu(SCN)2P and X=I 3 p respectively. For the a orbit in X=Cu(SCM)e a value of 3.3 ± 0.1 
has been obtained.ETEj Under a pressure of 20 kbar this effective mass decreases to 1.5m e £3 The rapid decrease 
in the effective mass with increasing pressure was correlated with the decrease in the superconducting transition 
temperature. A cyclotron resonance measurement (which should measure the bare mass and not the effective massEj) 
on X=Cu(SCN) 2 gives a mass of 1.2m e .El 

The measured effective masses are two to four times larger than predicted by Hiickel band-structure calculations 
that ignore the interactions between the electrons. Calculations using the local-density approximation (LDA) for 
X=Cu[N(CN) 2 ]Br give a band width a factor of two smaller than Hiickel.Ea 

The fact that the effective mass, specific heat, and magnetic susceptibility are significantly larger than predicted 
by band structure calculations suggest that the band width is narrowed by many-body effects. 

D. Superconducting phase 

There is increasing evidence that, like in the cupratesj^l the pairing of electrons in the superconducting state 
involves a different symmetry state than in conventional metals. Below I briefly review the evidence, pointing out 
where controversy exists. 

Nuclear magnetic resonance. The NMR relaxation rate 1/T\ does not have a Hebel-Slichter peak below T c . At low 
temperatures, 1/Ti ~ T 3 ^3 ! Ej instead of the exponentially activated behavior found in conventional superconductors. 
This power law behavior suggests there are nodes in the gap. 

Specific heat. At low temperatures the electronic specific heat of X=Cu[N(CN) 2 ]Br has a power law dependence on 
temperature, in contrast to the exponentially activated dependence of conventional superconductors. In a magnetic, 
field B perpendicular to the layers the electronic specific heat is linear in temperature with coefficient 7 ~ B 1 / 2 ^ 
This is the behavior predicted for a clean two-dimensional superconductor with nodes in the gap.Eil In contrast, in a 
conventional superconductor 7 ~ B. 

Magnetic penetration depth. There is controversy about its temperature dependence, dc magnetisation measure- 
ments on X=Cu£SiCN) 2 have been fitted to the exponentially activated temperature dependence expected for an s-wave 
superconductor .Ell In contrast, muon spin rotation (/+SR) measurements on X=Cu(SCN) 2 and X=Cu[N(CN) 2 ]Br show 
a linear temperature dependence with temperature.H Other /iSR measurements are consistent with an s-wave order 
parameter .Ej Surface impedance measurements on on X=Cu(SCNia, and X=Cu[N(CN) 2 ]Br give a penetration depth 
whose temperature dependence is consistent with s-wave pairingEj Furthermore, the real part of the conductivity 
versus temperature has a 'coherence' peak just below T c , as in conventional superconductors. 

Upper critical field H c2 . The H c2 versus temperature curve for X=Cu(SCN) 2 shows upper curvature near T c . The 
large slope near T c suggests that at low temperatures H c2 could exceed the Pauli limit (the field at which BCS theory 
predicts singlet Cooper pairs should be broken) .E30 

Vortex dynamics. Vibrating reed studies of X=Cu[N(CN) 2 ]Br in magnetic fields less than one tesla suggest anoma- 
lous vortex dynamics and the possibility of unconventional superconductivity.!! 2 ] 

II. THEORY 

I now argue that a theoretical model introduced by Kino and Fukuyamai captures enough of the essential physics 
of these materials to describe them, at least at the semi-quantitative level. Insight into the model is obtained by 
considering its various limits and what is known about similar models. 

A. Band structure 

The arrangement of the BEDT-TTF molecules in the k crystal structure and the dominant intermolecular hopping 
integrals are shown in Figure ||. Table I shows values for these integrals found from various quantum chemistry 
calculations for different anions. The intradimer hopping t^ is more than two times larger than the interdimer 
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integrals, t p and t q along the c + a direction, and tb 2 along the c direction. The dimer bonding and antibonding 
orbitals are split by approximately 2t;, 1 and so the mixing between bonding and antibonding orbitals can be neglected 
and the interdimer hoppings can be treated as a perturbation. Focusing on antibonding orbitals, each dimeric site 
has two nearest neighbors along the c direction (interaction t\ = if, 2 /2), and four next-nearest neighbors along the 
c± a directions (interaction ti = (t p + t q )/2). The dimer model is shown in Fig. ^ (compare Figures 1 and 6 in Ref. 

and Figure 13 in Ref. ) Note that each lattice site has hopping to six neighbours and so the lattice has the same 
co-ordination as the triangular lattice. Since each dimer has three electrons (or one hole), the upper (conduction) 
band is half-filled. 

If the interactions between electrons are neglected (which it will be argued below is not justified) then the model is 
just a tight-binding model and the dispersion relation for the half-filled band is 

E(k a , k c ) = 2t x cos (k c c) + 4t 2 cos (^y^J cos {^y^j (5) 

This band structure reproduces the main features of bands calculated using the extended Hiickel approximationE10. 
Bands calculated using the local-density approximation (LDA) for X=Cu(NCS)2 and X=Cu[N(CN)2]Br are also 
similar apart from an overall bandwidth narrowing.H The band structure (o) has been used to model the results of 
Shubnikov - de Haas oscillations measurementsE3 on X=Cu(NCS)2 and reflectance measurements^ on X=±3. 



B. Estimates of the Coulomb repulsion 



To test whether these materials can be described by a tight-binding model which neglects electron-electron inter- 
actions the effective Coulomb repulsion betweenJtwix.electrons in the anti-bonding orbital on a BEDT-TTF dimer is 
now estimated. Quantum chemistry calculationsLTLD estimate that the Coulomb repulsion between two electrons on 
a single BEDT-TTF molecule Uq is about 4-5 eV. These calculations neglect the screening effects present in the solid 
state and so the actual Uo may actually be smaller. Uo has been estimated from optical experiments to be 0.7 eV in 
(BEDT-TTF)HgBr3L3 and in a'-(BEDT-TTF) 2 Ag(CN) 2 to be 1.3 eV.El pp 

If Eo(N) denotes the ground state energy of N electrons on a dimer thenBD E (2) = \{U - (t/ 2 + WtlJ 1 / 2 ), 
E (3) = Uo — t, and E (A) = 2U . The effective Coulomb repulsion U between two electrons on a dimer is then 

U = E a (A) + E (2) - 2E (3) = 2t bl + \(U - (U* + 16^)1/2) (6) 

If Uq >> 4tf >1 the effective repulsion between two holes on a dimer is approximately 2t bl (See Figure ||). However, if 
Uq ~ ~ l eV then U ~ 0.3 eV. The optical conductivity data shown in Fig. g is consistent with this estimate. 
This value of U is comparable to the band width of the tight binding model and so suggests that correlation effects 
are important. 



C. A minimal model: the dimer Hubbard model 



Based on the above discussion wa-suggest that a minimal theoretical model to describe the k-(BEDT-TTF)2X 
crystals is the following HamiltonianjZl 

H = -h J2 ( c U c i,<r + h.c.) - t 2 (c| jCn# + />,c.) + P^(n lt - i)^ - ^) + /i^m,, (7) 

<ij>,cr <in>,o" i i,cr 

where c\ a creates a hole at dimer site i with spin projection a, Ui a is the hole number operator. The sum (ij) runs 
over pairs of nearest neighbor lattice sites in the horizontal direction, and the sum (in) runs over pairs of lattice sites 
along the diagonals. U is the on-site Coulombic repulsion on a single dimer (estimated in the previous section), and 
fi is the chemical potential. At half-filling \i may be non-zero due to the absence of particle-hole symmetry if both t\ 
and ti are non-zero. 

If there are strong correlations on each dimer then the amplitude for hopping between neighouring dimers, t\ and 
<2, are not simply given by t\ — i& 2 /2 and ti = (t p + t q )/2. The correct amplitude involves the overlap of a one hole 
state with the anti-bonding orbital creation operator acting on a two-hole state£il A straight-forward calculation gives 

— = - h - = -L (cos 9- sin9) (8) 
t b2 (t p + t q ) 2V2 
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where 




tan0=-£L- _1L +1 (9) 



If Uq >> 4ti, 1 then correlations on the dimer reduce the inter-dimer hopping by a factor of l/y/2. 

Various limits of the model (Q) are now considered leading to the speculative phase diagram shown in Fig. ^|. 

An important property of this model is that at half-filling a non-zero U is required for an insulating phase. A finite 
U is required because if t\ and t% are-both non-zero the nesting of the Fermi surface is not perfect. In the Hartree-Fock 
approximation Kino and Fukuyamafl found that for 2t\ = the critical value is U c /t2 — 3. 

If ti = t2 = t the model becomes the Hubbard model on a triangular lattice, after a rescaling of the lattice constants. 
Although there have been a number of studies of this model at half filling, no clear consensus of the phase diagram 
has been reached except that the ground state should be a paramagnetic metal for small U/t and an insulator for 
large U ft. i-The model has been studiepLby techniques including the Hartree-Fock approximation (which predicts 
U c /t ~ 5),&2l the slave boson technique^ and exact diagonalization of small clusters.^ The slave boson technique 
predicts that between the paramagnetic metal and insulating, antiferromagnetic phase that there is semimetallic linear 
commensurate spin-density-wave phase for 6.9 < U/t < 7.8.E3 In order to allow for the slightly different band structure 
associated with the lattice considered here these critical values of U must be divided by -\/3- Both the Hartree-Fock 
and slave-boson (mean-field) approximations tend to underestimate the effect of fluctuations and so the actual critical 
U may be larger than the above values. 

The estimate of U in the previous section and the parameter values given in Table I for various n— (BEDT-TTF^X 
crystals suggest that U / (ti + t^) ~ 3 — 6 and so these materials may be close to the metal-insulator transition. 



D. Strong coupling limit for the insulating phase 

If U 3> ti, t% then the ground state is an insulator and a standard strong-coupling expansion!! for the Hamiltonian 
implies that the spin degrees of freedom are described by the spin-^ Heisenberg model 

H = Jx J2 Si ■ Sj + Ji Si ■ S n (10) 

<ij> (in) 

where Si denotes a spin operator on site i, J\ = t\/U and J2 = t%/U. Again, the sum (ij) runs over pairs of nearest 
neighbor lattice sites in the horizontal direction, and the sum (in) runs over pairs of lattice sites along the diagonals. 
Ji and J2 are competing interactions leading to magnetic frustration. The parameter values in Table I suggest that 
J1/J2 ~ 0-3 — 1 and so magnetic frustration will play an important role in these materials. This model has sorpe 
similarites to a model considered for the layered titanate Na2Ti2Sb20 and with possible spin liquid ground states .EH 
Insight can be gained by considering various limits of this model. If Ji = then the model reduces to the Heisenbepg. 
model on a square lattice. At zero temperature there will be long range Neel order with a magnetic moment of 0.6/ib.l£I 
If J\ is non-zero but small it will introduce a small amount of magnetic frustration which will reduce the magnitude 
of the magnetic moment in the Neel state. The spin structure is showia-in Fig. ^ the same as that proposed for the 
antiferromagnetic phase of deuterated X=Cu[N(CN)2]Br based on nmrlia. As the ratio J\/ J% varies the wave vector 
associated with the antiferromagnetic order may also vary. There is the possibility of commensurate-incommensurate 
transitions. 

If Ji = J2 then the model reduces to the Heisenberg model on a triangular lattice. There has been some controversy 
about the ground state of this model. AndersonE3 originally suggested that the ground state was a "spin liquid" with 
no long-range magnetic order. However, recent numerical work suggests that there is long-range order but the quantum 
fluctuations are so large due tOr-the magnetic frustration that the magnetic moment may be an order of magnitude 
smaller than the classical valued 

If J2 = then— the model reduces to a set of decoupled Heisenberg antiferromagnetic chains which do not have 
long-range order £3 If J2 is non-zero but small the chains are weakly coupled. The case of only two chains corresponds 
to the "zig-zag" spin chain which is equivalent to a single chain with nearest-neigbour and next-nearest neighbour 
exchange, J 2 and Ji, respectively. (Unfortunately, the definition of ^v-4nd J 2 is the opposite of that normally used). 
This spin chain has been extensively studied and is well understood. E2l In the limit of interest here, J2 <§C Ji, there 
is a gap in the spectrum A ~ exp(— const. J1/J2) and there is long-range dimer order. We speculate that this "spin 
gap" is still present in the many chain limit. 
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Based on the above discussion a schematic phase diagram for the model (|7|) at half filling and at zero temperature 
has been drawn in Fig. ^. The possibility of superconductivity near the metal-insulator transition is suggested only 
on the basis of comparison with the experimental phase diagram. 

E. The effect of pressure 

A very important challenge is to understand how varying the pressure changes the ground state of the system. 
The experimental phase diagram shown in Fig. ^ certainly suggests that pressure reduces the electronic correlations. 
In most strongly correlated systems pressure increases the band width by more than the Coulomb repulsion and so 
reduces the effect of correlations. The effect of pressure in this system is more subtle because the Coulomb repulsion 
U is related to the hopping integral (see equation @) which will increase with pressure and the ratio ti/t 2 may 
vary with pressure and could drive the metal-insulator transition. Rahal et al.Ej measured the crystal structure of 
X=Cu(NCS)2 at room temperature and pressures of 1 bar and 7.5 kbar. They then performed Hiickel calculations 
of the hopping integrals for the corresponding crystal structures. The results, given in Table I, suggest that the 
amount of magnetic frustration decreases as the pressure increases. In contrast to the work of Rahal et al. a more 
heuristic approach was taken in an earlier study by Campos et al.E3 In order to fit experimental data for the pressure 
dependence of the Fermi surface area deduced from magneto-oscillation frequencies, they made assumptions about 
how the BEDT-TTF molecules move relative to one another as the pressure varies. They found an increase in ti/t 2 
with pressure. 

Assuming that the Coulomb repulsion on a single BEDT-TTF molecule Uq does not vary significantly with pressure 
then equation (|^) implies 

A Jh ^ (11) 

If Uo >> 4t6j then U will increase with pressure at about the same rate as which is similar to the rate at which 
t\ and t2 increase. This could mean that the correlations change little with pressure. However, if Uq ~ Atb t then the 
second term in ( pLl| ) reduces the rate at which U increases and so correlation effects can decrease with pressure. The 
uncertainty in the variation of the parameters with pressure is sufficiently great that more work is required before 
definitive answers can be obtained about this important question of how the ground state changes with pressure. 

F. Dynamical mean-field theory for the Hubbard model 

Some insight into the physics of the organics, particularly the existence of a low-energy scale in the metallic phase, 
can be gained by briefly reviewing progress that has been made in the past few years in understanding the Mott- 
Hubbard metal-insulator transition.E3 This has been done by studying the Hubbard model in a dynamical mean- 
field approximation which becomes exact in the limit of either large lattice connectivity or spatial dimensionality. 
This approximation maps the problem onto a single impurity problem that must be solved self-consistently. While 
time-dependent fluctuations are captured by this approximation spatially-dependent fluctuations are not captured. 
Some important physics that emerges from these studies is that there is an important low-energy scale To which is 
much smaller than the half-bandwidth D and the Coulomb repulsion U. This energy scale is the analogue of the 
Kondo temperature for the impurity problem. If there is significant magnetic frustration the Neel temperature for 
antiferromagnetic ordering is suppressed. Then at half filling and at fixed temperature T < T* ~ 0.051? as U is 
increased there is a first-order phase transition from a paramagnetic metallic phase to a paramagnetic insulating 
phase at U = U c ~ 3D. 

In the metallic phase the density of states p(u) contains peaks at energies ui = —U/2 and +U /2 which correspond to 
the lower and upper Hubbard bands, respectively and involve incoherent states. These peaks are broad and have width 
of order D. At temperatures less than To a new quasi-particle peak with width of order To forms at u> — (Figure 39 
in Reference |93|). The spectral weight of the peak vanishes as the metal-insulator transition is approached. The quasi- 
particle band involves coherent states that form a Fermi liquid. This means that the self energy S(o>) is independent 
of momentum has the energy dependence, ReS(o;) = U/2 + (1 — 1/Z)u> + 0(u> 3 ) and ImS(w) = —Bui 2 + 0{u A ). Z is 
the quasiparticle residue and is ceJated to the effective mass m* by m*/m e = 1/Z. As the metal-insulator transition is 
approached Tq — > and Z — > 0.0 The consequences of the above picture for the temperature dependence of different 
physical quantities for metals in close proximity to the metal-insulator transition are now briefly described. 



dU _ dt bl 
dP~ ~d~P 
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Optical conductivity^^ For temperatures T larger than To the optical conductivity is small near zero-frequency 
and has a broad peak at u> = U, corresponding to transitions between the lower and upper Hubbard bands. For T < To 
a Drude-like peak forms at zero- frequency and there is a new-shoulder at to ~ U/2. The latter peak corresponds to 
transitions from the quasi-particle band to the upper Hubbard band and from the lower Hubbard band to the quasi- 
particle band. The total spectral weight increases as the temperature decreases ..This picture has been used to explain 
the temperature dependence of the optical conductivity of vanadium oxides .EjE3o If this picture is used to explain 
the data shown in Fig. then the Coul omb repulsion U is estimated to be about 250 to 350 meV. This is roughly 



consistent with the estimates in Section [I B 



Resistivity. The temperature dependence of the dc resistivity p{T) is not monotonic. It has a peak at a temper ature 
of the order of To, and decreases rapidly below the peak. (See Figure 9 in Reference |9?], Figure 2 in Reference 100 , 
and the inset of Figure 2 in Reference |95| ). 

Mean-free path. Since coherent quasiparticles only exist for T < To it is not required that the mean-free path be 
larger than a lattice constant at high temperatures. Thus, it is possible for the resistivity to increase with temperature 
and be larger than the Mott minimum resistivity. 

Fermi liquid regime. The properties of the low temperature range (T < To) in which the metal has Fermi liquid 
properties is now considered. i_. 

Specific heat. In the Fermi liquid regime the linear coefficient iso 

7 = ^ (12) 
' 3ZD K ' 

where Z = 0.9(1 — U/U c ) is the quasi-particle weight. The effective Fermi energy is e* F = ZD. 
Effective mass. This is related to the quasiparticle weight Z by 

m* 1 
m e Z 

and diverges at the transition, consistent with the Brinkman-Rice scenaricll for the metal-insulator transition. 

Resistivity. At low temperatures p(T) ~ AT 2 , characteristic of a Fermi liquid. Near the Mott-Hubbard transition 
the coefficient A is related to the specific heat coefficient^ by A/-/ 2 — (2.3a) x 10~ 12 ohm cm (mol/mJ) 2 where a is 



the lattice constant in A of a three-dimensional system.E 

Magnetic susceptibilities. The uniform susceptibility (in units in which gps — 2) is 



*fe = °> = o354 + J (14) 

where J = D 2 /2U is the magnetic exchange energy. Note that this does not diverge at the transition. In contrast, 
the local susceptibility 

E. 4.7 . , 

x{q) = — (15) 



The local susceptibility can be related to the spin-echo decay rate 1/T 2 g 
The Wilson ratio. This is only universal for the local susceptibility. 
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/ free 

R= Xloc / X J oc =2.8 (16) 

where Xioc^ an( l 7^ ree are the values in the absence of interactions. 

The nmr relaxation rate. The spin relaxation rate 1/Ti depends linearly on temperature 

J- = h m ^M = i^ (17) 

TlT c^O LO (4) 2 v ' 

This could be combined with ( |l5|) to give a generalised Korringa-type relation. However, this is not directly relevant 
to nmr since the Knight shift K is proportional to the uniform susceptibility, not the local susceptibility. The Korringa 
ratio l/(TiTK 2 ) will diverge at the transition. 

I am unaware of calculations of the complete temperature dependence of 1/(T]T), the thermopowes, and Hall 



resistance. However, these quantities have been calculated for a doped Mott insulator. It was foundliHil that 1/Ti 
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versus temperature had a maximum at a temperature of the order of Tq. The thermopower versus temperature curvet 101 
has a maximum at a tenrnefafcme of the order of Tq at which its magnitude is of order 40 /iV/K. The Hallr*esistance 



is temperature dependent j!2d'E23 The Hall angle has an approximately quadratic temperature dependence.t°2L 

Luttinger's theorem is obeyed: the area of the Fermi surface is the same as in the absence of interactions.EU This 
means that even though there are strong electron-electron interactions, Hiickel band structure calculations (which 
ignore the interactions) may still be able to predict the Fermi surface area correctly. It has often been found for 
the organics, that Hiickel calculations can produce Fermi surface areas in agreement with the frequency of magneto- 
oscillations.B This success should not be used to argue that the electron-electron interactions are weak. 

In conclusion, many of the above properties are similar to those of the metallic phase of the k-(BEDT-TTF)2X 
salts, particularly the presence of a low-energy scale below which Fermi liquid behavior is observed. The dynamical 
mean-field theory results for temperature and frequency dependences for the conductivity appear to be very similar 
to experimental results. Two notable differences exist between the dynamical mean-field theory results and properties 
of the k-(BEDT-TTF)2X salts. First, the calculated value of the ratio A/j 2 is 5 to 200 times smaller than observed. 
Second, the nmr relaxation rate at low temperatures is not linear in temperature, as predicted by ([l7|). It remains 
to be seen if these discrepancies occur because for these quantities the two-dimensionality of the organics becomes 
important. 



III. FUTURE DIRECTIONS 



In conclusion, the experimental properties of this class of organic superconductors have been reviewed with an 
emphasis on similarities to the cuprates and how the metallic phase has properties quite distinct from a conventional 
metal. Quantum chemistry calculations can be used to justify a minimal microscopic model: a Hubbard model on 
an anistropic triangular lattice. Appropriate parameters for this model suggest that their is substantial magnetic 
frustration and the system is close to a metal-insulator transition. Consequently, recent results for the Mott-Hubbard 
transition suggest that this model can explain the unconventional metallic properties. However, there are many details 
to be worked out and so below some specific suggestions are made for future theoretical and experimental work. 

Theory. 

Based on the above discussion Figure ^ shows a speculative phase diagram for the Hubbard model (^) at half-filling 
and zero temperature. The challenge is to determine the actual phase diagram using the arsenal of theoretical tech- 



niques that have been developed over the past decade to attack Hubbard models far t, 
that could be used include exact diagonalizatiorEHj ajid, quantum Monte Casio li23'E££ 



le cuprates. Numerical techniqes 
Analytic techniqes that could be 



used include the fluctuation-exchapge approximation!^ and slave bosons.E£3 Insights might be gained from using idi 
from gauge theory of fluctuations!!!] and the SO(5) unified theory of superconductivity and antiferromagnetism.E^ 9 ] 
The most important question is whether this Hubbard model can produce superconductivity. This is controversial for 
the square lattice model away from half- filling J!23 Other questions of interest include: Does the metallic phase have 
non-Fermi liquid properties ? Ar e spin liquid phases possible? 

It was shown in Section II F that dynamical mean-field theory can provide important insights into the organics. 
However, the connection with the model Hamiltonian (|t]) needs to be made more precise by working out the appropriate 
lattice of large connectivity or dimensionality to study. Possibilities include the infinite connectivity limit of the 
Bethe lattice with nearest neighbour and next-nearest neighbour hopping and the two-sublattice frustrated model.Ej 
In both these models the degree of magnetic frustration can be varied. Then the complete temperature dependence 
of resistivity, Hall resistance, thermopower, magnetic susceptibility, and nmr relaxation rate should be calculated and 
compared to the organics. 

Finally, as discussed in Section [ll] a better understanding is needed of how the model parameters vary with pressure. 
Experiment. 

The organics have a number of advantages over the cuprates for experimental study. First, since pressure and not 
doping is used to tune through the metal-insulator transition (compare Figure ||) only a single sample is needed to study 
the transition. Second, the lower superconducting transition temperature and high purity samples mean that magneto- 
oscillations provide a powerful probe of the Fermi surface. On the other hand, the single crystals are sufficieatly small 
that neutron scattering is difficult and the first studies on the organics have only recently been performed." 1 . 1 ! Further 
studies could probe the magnetic fluctuations. X=Cu[N(CN)2]Br lies close to the metal-insulator transtionc3 and so 
is the ideal candidate for further study. More systematic experimental studies of the metallic state and its variation 
with pressure are needed. Extensive transport and thermodynamic measurements should be made on the same single 
crystal__.Of particular interest is how much of the low-temperature data can be explained within a Brinkman-Rice 
picturd_j of a Fermi liquid with an effective mass that diverges at the metal-insulator transition. 
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By applying uniaxial stress along the c-axis (in the conducting plane) it shotuld .be. possible to vary the ratio of t\ 
to t% and thus vary the magnetic frustration. A few uniaxial stress experimentslll5tij have been performed. One did 



find that for stress in the plane the resistivity peak and superconducting transition temperature were enhanced.t!2 
More extensive measurements are needed. 

X=Cu2(CN)3 is particularly interesting because-in the insulating phase at ambient pressure there is no evidence 
of antiferromagnet ordering at low temperatures £3 This should be checked with nuclear magnetic resonance. The 
magnetic susceptibility decreases rapidly below 10 K. Is this due to a spin liquid ground state? The parameter values 
in Table I suggest that this material has the largest amount of magnetic frustration of the listed materials. 
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Note added in proof. Since submission of this article preprints have appeared which consider the possibility of super- 
conductivity in the Hubbard model on the anisotro pic triangular lattice . Calculations at the level of the random-phase 
approximation (RPA) [M. Vot.ja a nd E. Dagotto, cond-mat/9807168[ and the fl uctuation-exchange approxi mation 
fj. Schmalian , |cond-mat/9807042fc H. Kino and H. Kontani, |cond-mat/9807147t H. Kondo and T. Moriya, |cond ' 



mat /9807322| ] suggest that at the boundary of the antiferromagnetic phase the model exhibits superconductivity 
mediated by spin fluctuations. This result is consistent with the speculative phase diagram shown in Figure 6 of this 
article. As the ratio t\/t2 increases the wavevector associated with the antiferromagnetic spin fluctuations changes 
from (-7T, 7r) and the RPA calculations suggest that the superconductivity changes from d-wave singlet to s-wave triplet 
in the odd-frequency channel. 
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FIG. 1. Schematic phase diagram of the k— (BEDT-TTF)2X family of organic conductors 
antiferromagnetic phase, PMI an paramagnetic insulating phase, and SC a superconducting phase 



' ll6 l AFI denotes an insulating 
The phase transition 



between the paramagnetic phases is first order and denoted by the dashed line which ends at a second order critical point. 
The arrows denote the location of at ambient pressure of materials with different anions X. As the pressure increases the 
unconventional properties of the metallic phase disappear. The above phase diagram is qualitatively similar to that of the 
cuprate superconductors with doping playing the role of pressure. 
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FIG. 2. Non-Drude behavior of the frequency-dependent conductivity <r{u>) of X=Cu[N(CN)2]Br. The data is from Reference 
^l] and shows the strong temperature dependence of the low- frequency conductivity. At 25 K and below (t(lu) has a Drude-type 
peak at zero frequency (see the dashed curve). In contrast, this peak is completely absent at 50 K and above (see the solid 
curve). The broad peak around 300 meV can be identified with transitions between the lower and upper Hubbard bands. The 
very sharp features are due to infra-red active phonons. The data shown is for the electric field parallel] to the a-axis. Similar 
data is obtained for the electric field parallel to the c-axis except the broad peak is around 400 meV.Eil (Both the a-axis and 
the c-axis are in the plane of the layers). 
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FIG. 3. Geometrical arrangement of the BEDT-TTF molecules in the conducting layers of the k phase. Each line represents 
a BEDT-TTF molecule. The dominant intermolecular hopping integrals are shown. Some calculated values for these integrals 
are given in Table 1. 




FIG. 4. Triangular lattice for the dimer model of k— (BEDT-TTF) 2 X. Each lattice site is denoted by a circle and represents 
the anti-bonding orbital on a dimer, a pair of BEDT-TTF molecules. Each dashed oval represents a single BEDT-TTF molecule. 
There is no hopping in the vertical direction and there is hopping along the diagonals. Competition between the two types 
of hopping results in magnetic frustration and a Mott-Hubbard metal-insulator transition at a critical value of the Coulomb 
repulsion. 
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FIG. 5. Effective Coulomb repulsion between two electrons in the anti-bonding orbital of a dimer in the strongly correlated 
limit, (a) The ground state of a pair of three-quarter filled dimers. If tt 1 is the intradimer hopping integral then 2tt 1 is the 
energy splitting of the bonding and anti-bonding orbital, (b) If the Coulomb repulsion Uo on a single molecule is much larger 
than tb 1 then this denotes the configuration of the lowest energy configuration involving charge transfer between dimers. The 
energy difference between (a) and (b) is 2tt 1 and so this is the effective Coulomb repulsion on a dimer for Uo ^> tb 1 . 
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FIG. 6. Speculative phase diagram for the Hubbard model on an anisotropic triangular lattice (|7|). Shaded and unshaded 
areas represent insulating and metallic regions respectively. The solid line with an arrow shows how the ground state of a 
particular k— (BEDT-TTF)2X crystal might change as the pressure is increased. Vertical lines on the left side, center, and 
right side correspond to the Hubbard model on a square lattice, triangular lattice and decoupled chains, respectively. 
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FIG. 7. Collinear spin structure of the antiferromagnetic ground state for the dimer model ([7|) with t\ <C ti <C (7. Each circle 
represents a dimer of BEDTptTF molecules. This is the same structure as that proposed for deuterated X=Cu[N(CN)2]Br 
based on nmr measurements. E-3 
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FIG. 8. A simple unified picture of the phase diagram of the organics and cuprates at low-temperatures. The organic super- 
conductors k-(BEDT-TTF)2X are always at half- filling and one tunes through the insulator-superconductor-metal transition 
by applying pressure. In the cuprates this transition is passed through by varying the band filling. 
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TABLE I. Values of the hopping integrals (in meV) between the BEDT-TTF molecules in the k-(BEDT-TTF) 2 X crystal 
structure for different anions X. All values were based on quantum chemistry calculations using the Hiickel approximation, 
except the first line which involved an ab initio calculation. All values are for the crystal structure at ambient pressure, unless 
denoted otherwise. The parameters in the Hubbard model ([?]) are given by t\ ~ tb 2 /2, fa ~ (tp + t q )/2, and U ~ 2tf >1 (See 
equations (|^) and (|^)). The fact that ti ~ tu implies that there is significant magnetic frustration. Note that for the same 
anion there is some variation between the values of the hopping integrals calculated by different researchers. 



Anion X 




tb 2 


tp 


t q 




Ref. 


Cu[N(CN) 2 ]Br 


272 


85 


130 


40 


0.50 


mm 


Cu[N(CN) 2 ]Br 


224 


71 


94 


40 


0.53 


■i 


Cu[N(CN) 2 ]Br 


244 


92 


101 


34 


0.68 


■e 


Cu(CN) 3 


224 


115 


80 


29 


1.10 


■e 


Is 


247 


88 


119 


33 


0.58 


l E 


Cu(SCN) 2 


230 


113 


99 


33 


0.85 


ll 


Cu(SCN) 2 


136 


50 


41 


27 


0.74 


mm 


Cu(SCN) 2 (7.5 kbar) 


162 


54 


58 


30 


0.61 


Mm 


Cu(SCN) 2 


244 


26 


22 


20 


0.62 




Cu(SCN) 2 (20 kbar) 


324 


45 


27 


31 


0.88 
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